Kondo effect and bistability in a double-quantum-dot 
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We study theoretically the out-of-equilibrium transport 
properties of a double quantum dot system in the Kondo 
regime. We model the system by means of a two-impurity 
Anderson Hamiltonian. The transport properties are char- 
acterized by Kondo effect properties, however, superimposed 
them, the system possesses novel non-linear bistability behav- 



Recently experiments on quantum dots (QDs) at tem- 
peratures (T) below the Kondo temperature (T^) have 
shown that new physiGa|-|emerges when their transport 
properties are studied-BB ' These experiments confirm 
that many of the phenomena that characterize strongly 
correlated metals and insulators, as it is the case of the 
Kondo effect, are present in QDs. The advantage of 
studying the quantum-coherent many-body Kondo state 
in QDs, in comparison with natural compounds, con- 
sists in the possibility of continuous tuning the relevant 
parameters governing the properties of this state, in an 
equilibrium and out-of-equilibrium situation. In particu- 
lar, the problem of electrons tunnelling through double- 
quantum-dots (DQDs) in the Kpndo regime has received 
much attention in recent years OH The DQD is the sim- 
plest system where it is possible to study the compe- 
tition between the inter-dot antiferromagnetic spin-spin 
correlation and the dot-conduction spin-spin correlation 
present in its ground state. The type of coupling between 
the QDs determines the character of the electronic states 
and the transport properties of the artificial molecule. In 
the tunnelling regime, the electronics states are extended 
across the entire system and form a coherent state based 
on the bonding or anti-bonding levels of the QDs. In 
this context, recently Aguado and LangrethQ studied the 
out-of-equilibrium transport properties of a DQD in the 
Kondo regime. They found that for inter-dot-coupling 
greater than the level broadening, there is a critical volt- 
age above which the coherent configuration is unstable. 
This instability is reflected as a drastic drop of the cur- 
rent leading to a singular region of negative differential 
conductance. This behavior resembles the I-V charac- 
teristic of a double-|barrier structure in the accumulation 
of charge regime. tj This system, due to non-linearities 
introduced by the Coulomb interactions, has a bistable 
behavior, characterized by two solutions for the current. 
In this work we report the existence of a similar bistable 
behavior in an out-of-equilibrium double quantum-dot in 
the Kondo regime. 



We focus our study on a DQD in the Kondo regime 
driven to an out-of-equilibrium state by means of a dc 
voltage bias. In recent years the very small width of 
fabricated semiconductor leads used to connect the QDs 
permits to separate each lateral quantum confine energy 
level. In this case the leads can be represented by a 
ID tight-binding Hamiltonian connecting the dots to two 
particle reservoirs characterized by Fermi levels /il and 
fiR respectively. It has been shown that the physics asso- 
ciated to a DQD connected to two leads can be readily 
understood in terms of two impurities Andgrson Hamil- 
tonian where the impurities are the QZJsOqlI The Hamil- 
tonian can be written as, 
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where the dots have been localized at sites and 1 of the 
lattice. The operator c\^ creates an electron in the site 
i with spin a, Si is the site energy, t is the hopping in 
the leads, Vj;,(ij) is the hopping between the left (right) 
leads and the left (right) QD, tc is the inter-dot coupfing 
tunnelling and U is the on-site Coulomb energy. In typ- 
ical experiments, the intra-dot Coulomb repulsion U is 
large compared with kT. In this case, taking the limit 
[/ ^ oo, it is possible to obtain an important simplifica- 
tion of the model because the double occupancy at the 
dots is eliminated from the Hilbert space. 

The Hamiltonian may be written in tp;raa of auxiliary 
slave boson operators plus constraints. EiJtJThe annihi- 
lation operator of an electron at dot a is decomposed as 
Cq(t = bl^faa, where a slave boson operator 6j, creates 
an empty state and a fermion operator faa- annihilates 
a single occupied state with spin a. The sites a = 0, 1 
can only be in one of the states 6]^ |0) or f^^ |0). The 
slaves boson operator acts preventing double occupancy 
of the site. When an electron creation operator acts on 
an occupied site, the boson part annihilates the state; 

4,. fL |o) = /l.fea fL |o) = 0. 

The exclusion of double occupancy at site a, imposes 
the condition, Qa = J2a flafacr + bl^ba = 1, for the 



(1) 



number of bosons and fermions at that site.ll3 This con- 
straint can be taken into account by adding a term, 
J2a ^a{Qa ^1)) to the Hamiltonian (Eq.l) with the La- 
grange multiphcrs A^that are calculated imposing the 
constraint condition. In order to solve the Hamiltonian 
we adopt the mean field approximation (MFA) where the 
Bose operators 6jj and ba are replaced byJie expecta- 
tion values, (ba) = baV^ = (ba) = 6j,'\/2,t2l neglecting 
its fluctuations. 

Adopting the slave-bosons MFA, the tight-binding 
two-impurity Anderson Hamiltonian describing a double 
quantum-dot connected to leads can be written as, 

i <ij=iO,l>cr Q=0,1,<T 

a a 

+ ^cE(/ol/i- + /iUo.) + E A.(6t.6„-1). (2) 



where £o_= £o + Ao, £i = £i + Ai, Vl = Vl^o, ^r = Vr&i, 
tc = tcbobi. 

On a tight-binding basis, the stationary state of energy 
Ek results to be 



\^ka) =J2"'t\(t>i. 



(3) 



where \4>ia) is a Wannier state localized at site i of spin 
a and the coefficients a'^^ obey the non-linear difference 
equations, 

Ekal, = Ejal^ + t{a]_^^, + a^^i,,) (j ^ -1,0, 1,2), 

A; k -vr k J. k 

^fe«0(l),cr = £0(1)00(1),^ + VL{R)a-l(2),a + ^ca^Qy^, (4) 



where 



J.o- 



is the amplitude of probability to find the 



electron in the site j in the state k with spin a. The 
mean field parameters (60, 61, Aq, Ai) are determined by 
minimizing the free energy of the systems. Taking the 
expectation value of the Hamiltonian, differentiating and 
using the Hellman-Feymann theoremJ13it is possible to 
find that they satisfy the set of four equations: 
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The sum over the wave vector k cover all the occupied 
states. The resulting equations are nonlinear because of 
the renormalization of the localized levels in the dots, the 



inter-dots coupling tunneling and the coupling tunneling 
between the QDs and the leads. Due to the presence of 
these non-linear terms, the set of equations (4) and (5) 
have to be solved self-consistently. 

In order to study the solutions of equations (4) and (5) 
we assume that the electrons are described by a plane 
wave incident from the far left with an amplitude / and 
a reflection amplitude R and at the far right by a trans- 
mission amplitude T. Taking this to be the solution we 
can write, 
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where Vj is the position of site j. The solution of equa- 
tions (4) can be obtained through an adequate iteration 
of it from right to left. For a given transmitted am- 
plitude, the associated reflected and incident amplitudes 
may be determined by matching the iterated function to 
the proper plane wave at the far left. The transmission 
probability \T/I\ , obtained from the iterative procedure, 
multiplied by the wave vector k, gives us the contribu- 
tion of this wave vector to the current. With the pur- 
pose of solving equations (4) an (5) we define a second 
pseudo-time-like iteration in the following way. Initially 
Eqs. (4) are solved, for a particular applied potential, 
V, ignoring the non-linear term and for energies up to 
the maximum Fermi energy. The coefficients thus ob- 
tained correspond to a solution for non-interacting elec- 
trons. They are used to construct the non-linear term 
for the next solution. The procedure is repeated up to 
the moment in which convergence is reached. This so- 
lution is taken as the starting point for the next cycle 
corresponding to another value of V. 

Once the amplitudes a^^ are known, the current is 
numerically obtained from. 
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We study a model which consists of two leads con- 
nected to quantum dots with fi^ = —V/2 and fj^n = V/2, 
t = 30 F, Vi = Fi? = Vb = 5.48 F, eq = £1 = -3.5 F 
(Kondo regime with Tk « IQ-^ F with F = nV^piQ)). 
The normalization of the wave function is taken so that 
each site of the leads can be populated by a maximum of 
two electrons. 

Fig 1 shows the I-V characteristics. The system has 
a bistable behavior, with two solutions for the current, 
within a range of values for the external potential that in- 
creases with tc ■ These solution are obtained as the volt- 
age is increased (solid line), or decreased (dashed line). 
When it is increased above a critical value Vc^, the ccb: 
herence between the dots is lost and the current drops .Q 
On the other hand, when it is reduced from above Vc^ , 
the Kondo resonances in each dot are pinned at their 
own chemical potential and it is necessary to reduce the 



voltage by an additional amount AV, below Vc^ , to re- 
store the coherence between the dots. The width of 
the bistability AV is proportional to tc and it is about 
Ay ^ lijLeV{ tc — 1.3T)~. that is within the resolution 
limits of the experimentsu. 

It is well known that the slave-boson MFA describes 
a system that, for a region of the parameters space or for 
T > Tk, possess a non-physical solution that disconnects 
the magnetic impurity from the rest of the system. For 
the problem we are studying of a DQD this solution co- 
exist with the physical onctj. This non-physical behavior 
is not associated to the bistability we are describing. To 
make this point clear in Fig.l we represent the spurious 
solution (dot line) , as a function of the external potential 
with J = for y < V/, / 7^ 0, ^^ > V* . 

For ic > r, the Kondo resonances are split into a 
bonding and antibonding peaks. They are located at, 

e± = (eo + ei) ± \J {e^ -eif + Atl}/2. In fig 2 we plot 
£± and £0(1) ^-s a function of V for tc = l.SF. For 

V = Q, the bonding and antibonding levels satisfies e± 
= itcand the site energies £o(i) = 0- When the bias V 
is increased, e± are kept almost constant and eo and ei 
decreases and increases respectively. For higher values of 

V (always V < Vc^) ^ £+ and e_ approach each other, 
?o goes up and ei down, until they almost coincide again 
at y = Vc^. For a voltage close to and above K^, e± 
and £0(1) converge towards their own chemical potential 
f^L(R)- When the external potential is decreased from 
above Vc^ , the £± are kept near their own chemical po- 
tential until a critical voltage Vc, is reached, where they 
recover abruptly the values obtained when the potential 
was increased from zero. Decreasing the voltage in the 
interval Vc^ < V < Vc^, the separation between the £± 
is less than V and reduces with V until V = Vc^ ■ Below 
Vc^ , the system returns to a coherent state and the e'± 
go again to the neighborhood of ±ic- 

Fig 3 shows the left and right DOS of the QDs for 
tc = l.Sr and a voltage inside the bistability region 
{V = 2.15Tfc). The figure shows a coherent solution with 
the dots strongly coupled (high current; solid line), with 
an equally shared split Kondo peak and another one with 
the dots weakly coupled (low current; dashed line), each 
having a DOS with a Kondo peak pinned at their corre- 
sponding Fermi energies. 

The nonlinear behavior described above has similari- 
ties with the well-known bistability present in 3D reso- 
nant tunneling double bamex-sjistems and multistabili- 
ties in doped superlattices£3'Ej~llj. In 3D tunneling dou- 
ble barrier systems, the Coulomb interaction pins the 
renormalized energy level at the well when the applied 
voltage V is augmented maintaining the system at res- 
onance. However, when V goes above a critical value 
the abrupt leakage of the charge accumulated in the well 
takes the system out of resonance. For doped super- 
lattices, at large carrier densities, the internal poten- 
tial profile produced by the applied voltage breaks up 



into dopiaiHS giving rise to a multistable behavior in the 
currentlljll3. In both systems, the non-linear deforma- 
tion of the potential profile due to the electron-electron 
interaction is responsible for the phenomenon. Although 
in our case the Kondo effect is a strong many body ef- 
fect of a different nature in comparison with the pre- 
vious systems, the bistability is as well a consequence 
of the non-linear potential profile deformations that the 
electrons suffer when the external potential is modified. 
Increasing the voltage (always V < Vc^ ), the interac- 
tion maintains the site energy of the right dot above the 
left dot level as shown in Fig 2. For these values of V, 
the electron wave function delocalizes and spreads on the 
double-quantum dot. In contrast just above T4^, this 
configuration becomes unstable, the coherence between 
the dots is lost, their interaction reduces abruptly as the 
external potential is augmented and the wave function 
localizes at each quantum dot. The tunneling coupling 
between the dots {tc) reduces abruptly at Vc, and con- 
sequently the current circulating along the system (Fig 
4). For the DQD as it was as well for the double barrier 
heterostructure and the superlattices, due to the nonlin- 
earities produced by the Coulomb interaction, this tran- 
sition occurs at different values of V, depending wether 
the potential is increased or decreased, giving rise to the 
bistability structure present in the I-V characteristics. 

In conclusion, we have reported that a double 
quantum-dot could have a bistable behavior in the Kondo 
regime at zero temperature. This evidence was obtained 
using a mean-field approach that neglects fiuctuations, 
appropriate for large spins as it is the first term of an 
1/A'^ expansion, where -/V-iis the spin dimension. Some 
experiments have shownn3 that it is possible to design 
dots with a configuration that, due to the exchange in- 
teraction, the parallel spin coupling following Hund's rule 
gives rise to a total dot spin St > 1/2. For this cases 
our results are more reliable as fluctuations are less im- 
portant. Although we believe that our results are robust 
against them, a study that goes beyond mean field is nec- 
essary to confirm our conclusions. As mentioned above 
the slave boson MFA have a non-physical solution and 
it is possible to speculate that the negative differential 
resistance and the bistable behavior of the current could 
be associated to it. However, we were able to show that 
the phenomena here studied are not associated to this ar- 
tifact of the mean field approximation as they correspond 
to different solutions of the problem. 

Although bistabilities are a well-known phenomenon in 
resonant tunneling, to the best of our knowledge neither 
its existence has been proposed nor has it been measured 
in a ID Kondo system. This phenomena is expected to 
occur due to the high nonlinearity of the problem. Other 
possible manifestations arising from nonlinearity in DQD 
are multistabilities, self-sustained oscillations and chaos 
phenomena, that are currently been investigated. 
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FIG. 4. a) Renormalized level broadening in the left and 
the right dots as function of the voltage, b) the renormalized 
interdot coupling tunneling versus voltage for tc — l.SF 
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FIG. 1. I - V curve for a) tc = l.OF, h)tc = I.IF, 
c)tc = 1.2r and d) tc = 1.3r 



FIG. 2. Level energies, a) e± vs V and b) eo(i) vs V for 

tc = i.3r 



FIG. 3. Local density of states for left and righ dot for 
tc = 1.3 and V = 2.15Tfc 
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